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CNl ' Abstract 
> 

' The rare decays B® — * D^n T and B° — * D^tt^ can occur only via annihilation 

type diagrams in the standard model. In this paper, we calculate branching ratios 

o ' 
in 
o 

Q_i| asymmetry, which may be measured in LHC-b experiment in future. 

1 Introduction 

• rH - 

X 

The rich data from two B factories make the study of B physics a very hot topic. A lot 
of study has been made, especially for the CP violation problem. The CKM angle (3 = 0i 
has already been measured pQ. However the other two angles are difficult to measure in B 
factories. The study of B s meson decay is needed for this purpose. Some work on B s decays 
have already been done [21 E] . 

In this work, we will explore four decay channels, namely B° s — > D^ir^ and — > D T 7i ± . 
There is only one kind of contribution for each of the decay mode, thus there is no direct 



of these decays in perturbative QCD approach ignoring soft final state interaction. 
From our calculation, we find that their branching ratios are at O(10~ e ) with large CP 
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CP violation for them. However there is still CP violation induced by mixing, although they 
are decays with charged final states (non CP eigenstates). They are quite complicated since 
altogether four are involved simultaneously [I]. 

From these decays, we find that the four quarks in final states are different from the ones 
in B° s meson. We call this mode pure annihilation type decay. In the factorization approach, 
this decay is described as B° s annihilating into vacuum and final states mesons produced from 
vacuum afterwards. They are rare decays. Up to now, only PQCD approach can calculate 
this kind of modes effectively. Using PQCD approach, we have calculated many of this kind 
of decays [SI El E], and some decays have been measured in B factory. Some information 
about PQCD in detail can be found in ref.[7]. 

In standard model language, for decay B s — > Dtt, a W boson exchange causes bs — > uc, 
and dd in final state are produced from a gluon. This is also called W exchange diagram. 
This gluon can be emitted by any one of quarks participating in the four quarks interaction. 
This is shown in FigureJU We consider the B® meson at rest for simplicity. In this frame, 
this gluon has 0(M B /2) momenta, that's to say, this is a hard gluon. We can perturbatively 
treat the process by six quarks interaction. 
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Figure 1: The picture of PQCD approach. 



In this work, we will give the PQCD calculation of these two decays in the next section, 
and discuss the numerical results in section [31 At last we conclude this study in section HI 
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2 Calculation 



The non-leptonic B® decays B® — > D + n and B® — > £) 7r + are rare decays. For decay 
S° — > D + n~ , the effective Hamiltonian at the scale lower than My/ is given j8] as: 

Wi = % [c^c^O) + c 2 (/i)o 2 (/i)] , (f ) 

where the four-quark operators are 

Ol = (fo)y_ A (cM)y_ A , 2 = {bu) V - A {cs) V -A , (2) 

with the definition (gig^v-A = 9i 7/^(1 — 75)0*2- Ci,2 are Wilson coefficients at renormalization 
scale jji,. For decay 5° — > D~tt + , the effective Hamiltonian read: 

where the four-quark operators are 



cimqm + c'mo'm , (3) 



0[ = (bs) v _ A (uc) v _ A , 0' 2 = (bc) v ^ A {us) v _ A . (4) 

In PQCD, the decay amplitude is expressed as [7j: 

Amplitude - [d%d A k 2 d 4 k 3 Tr [C(t)$ B o(ife 1 )$ J j(A; 2 )$ w (A;3)^(/: 1 , fc 2 , fc 3 , t)e" 5(i) ] . (5) 



In this Equation, C(t) is Wilson coefficient at scale t with leading order QCD correction. 
are light-cone wave functions, which describe the non-perturbative contributions. They can 
not be theoretically calculated directly. Fortunately, they are process independent. e~ s ® is 
called Sudakov factor, which comes from the resummation of soft and collinear divergence. 
This Sudakov factor suppress the soft contributions, which make the perturbative calculation 
of hard part reliable. By including the dependence of the wave functions and Sudakov 
form factor, this approach is free of endpoint singularity. Thus, the work left is calculating 
the perturbative hard part Hit). 

The structures of the meson wave functions are 

B° S (P): [f + msftsMz), (6) 
D(P): >y s [f + m D ]<l> D (x), (7) 
tt(P) : ^<s[P(j) A (x) + m (f)p(x) + m (j6 jh - l)4> T (x)] , (8) 



with m = m 2 /(m u + ma) = 1.4 GeV, characterizing the chiral breaking scale. And the 
light-like vectors are defined as n = (1, 0, Oy) and v = (0, 1, 0^). In above functions, 4>i are 
distribution amplitude wave functions. 

According to the effective Hamiltonian 01) , the diagrams contributing to B° s — > D + n~ 
are drawn in Figf2J Just as stated in Section HJ this decay has only annihilation diagrams. 
With the meson wave functions and Sudakov factors, the hard amplitude for factorizable 
annihilation diagrams in Figf5]^a) and (b) is, 

F a = —z-Mb \ dx 2 dx 3 / b 2 db 2 b 3 db 3 (f) D (x 2l b 2 ) x {(1 - x 3 ) 4>^{x 3 ) 
-J Jo Jo 1 

+ r (3 - 2x 3 )r n (f)P(x 3 ) - r(l - 2x 3 )r 7T (f)l(x 3 )}E f (t 1 a )h a (x 2 ,x 3 ,b 2 ,b 3 ) 

- {x 2 ^{x 3 ) + 2r (1 + x 2 )r^ (x 3 ) }E f (t 2 a )h a {l -x 3 ,l-x 2 , b 3 , 6 2 )] , (9) 

where r = m£>/M B o, r v = m Q /M B o, and the functions Ef containing Sudakov factors and 
Wilson coefficients of four quark operator, hard scale t]; 2 and virtual quark and gluon prop- 
agator h a are given in the appendix. The explicit form for the distribution amplitude (fi M 
of wave functions, are given in the next section. In above function, Xj denotes light (anti-) 
quark momentum fraction in meson. 

For the non-factorizable annihilation diagrams in Figj2]Jc) and (d), the results can read: 

2567T f 1 [°° 
M a = — ==M| / dx x dx 2 dx 3 \ b x db\b 2 db 2 4 , B o(xi,b 1 )(f) D (x 2 ,b 2 ) 
6\/lJS c Jo J 

x {x 2 (f>^(x 3 , b 2 ) + r (1 + x 2 - x 3 ) r n <f>%(x 3 , b 2 ) 
+ r (1 - x 2 - x 3 ) r n (f)l(x 3 , b 2 )}E m (t l m )h ( a l \x 1 , x 2 , x 3 , h, b 2 ) 
- { (1 - x 3 ) 4>^(x 3 , b 2 ) + r (3 + x 2 - x 3 ) r n (f)^(x 3 , b 2 ) 

+ r(x 2 -l + x 3 ) r w <f>l(x 3 , b 2 )}E m (t 2 JhW(x 1 ,x 2 , x 3 , h, 6 2 )] . (10) 

The total decay amplitude for _B° — > D + n~ is given as 

A = f B F a + M a . (11) 

The decay width is expressed as 

T(B° - = ^(1 - r 2 )\V: b V csA \ 2 . (12) 
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Figure 2: Leading order Feynman diagrams contributing to decay B® — > D + 7i 



As the case B° s — > D~n + , we also draw diagrams Fig J3] using Equation (jHJ). The amplitude 
for factorizable annihilation diagrams (a) and (b) results in —F a . The amplitude for the non- 
factorizable annihilation diagram results in 
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M 2 B \ dx\dx 2 dx 3 \ b\db\b 2 db 2 4>b°{xi, &i)0z?(^2, ^2) 



{(1 - x 3 ) <p^(x 3 , b 2 ) + r (x 2 + 1 - x 3 ) r w (f)^(x 3 , b 2 ) 
+ r(x 2 -l + x 3 )r w (frl(x 3 ,b 2 )}E m (t l m )h£ ) (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 
- {x 2 (p^(x 3 , b 2 ) + r(3 + x 2 - rr3)rv</£(x 3 , fe 2 ) 

+ r(l - x 2 - X3)r 7r ^(x3,6 2 )}-Em(^) /l i 2) ( a; i ? a; 2,a;3,&i,&2) • (13) 
Thus, the total decay amplitude A' and decay width V for B° s — ► D~7T + decay is given as 

A' = -/ b F + M;, (14) 
ZT7T+) = ^^(l-r 2 )IO^f. (15) 



r(5 s L 
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The decays widths for CP conjugated mode, B® — > D T 7r ± , are the same expressions as 
5° — > D ± 7r =F , with the conjugate of CKM matrix elements. 
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Figure 3: Leading order Feynman diagrams contributing to decay — > D ir 



3 Numerical Evaluation 



Considering SU(3) symmetry, we use the distribution amplitude of the meson similar to 
B meson: 



4>b°( x , b) = Nx 2 (l — x) 2 exp 



2a;? 2 V ; 



(16) 



which is adopted in ref. [SI QUI HI]- iV is a normalization factor, which can be get from 
normalized relation: 



Gfe (f) M (x,b = 0) = 
Z\fZl\ c 



For .D meson, the distribution amplitude is 

, , > 3 



f D x(l-x){l + a D (l-2x)}, 



(17) 



;is) 
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which is fitted from experiments [12] . The wave functions of the 7r meson have been derived 
in ref.pSJQI]: 

3/. 



y/2N L 



2v/2AT 



with the Gegenbauer polynomials, 



x(l - x) 1 + 0.44C 2 3/2 (2x - 1) + 0.25C4 /2 (2x - 1) 
1 + 0.43C 2 1/2 (2a; - 1) + 0.09C 4 1/2 (2a; - 1) 
1 - 2x) [I + 0.55(10x 2 - lOx + 1)] , 



(19) 
(20) 
(21) 



Cl l2 (t) = l(3t 2 - 1) , Cl /2 (t) = i(35t 4 - 30t 2 + 3) , 
Cl /2 {t) = - 1) , Ct /2 (t) = ^(21t 4 - 14t 2 + 1) . 



(22) 



The other input parameters are listed below [15J: 



= 230 MeV , u B = 0.5 GeV , f D = 240 MeV , C D = 0.8 ± 0.2 , f n = 132 MeV , 
m B = 5.37 GeV , m D = 1.87 GeV , m = 1.4 GeV , t b o = 1.46 x KT 12 s , 
\V cb \ = 0.043, = 0.22 , 1^1= 0.0036, |K S | = 0.974 . (23) 

With above parameters, we show the decay amplitudes calculated in Tabled The pre- 
dicted branching ratios are: 



Br(£° 
Br(B° s 



D + 7i- 

ZTtT 



13 x 10 



2.9 x 10 



-7 



-6 



(24) 
(25) 



From above results, we find that the branching ratios of decay B° s — > D + ir is smaller 
than that of decay B° s -> D-n + . The CKM element in decay B° s -> D+ti- is V* b V cs , but 
in B° s -> the CKM element is V* b V us . Although \V* h V us \ and \V* h V cs \ are both C(A 3 ) 

in Wolfenstein parametrization, the value of ^"A^" C3 | is equal to 0.37. The branching ratio of 



\v: b v us \ 

3° 

s 



Bg — * D~7i + is 3 times larger than that of B® — ► D + n~, which is mainly due to the CKM 
factor. 

In addition to the perturbative annihilation contributions, there are another pictures 
existing such as soft final states interaction [16] . In ref.[5], the results from the PQCD 
approach for B° — > DjK + is consistent with experiment well, which tells us the soft final 



Table 1: Decay amplitudes (10 3 GeV) with parameters eqs. (!9lfT31. 



B° s 


-> D + TT- 


B° s 




f B F a 


0.51 - 1.3z 


-f B F a 


-0.51 + 1.3i 


M a 


-16.1 - 19.1 i 


M' 

a 


-1.8 - 19.1 i 


A 


-15.6 - 20Ai 


A' 


-2.3 - 17.8 % 


Br 


8.3 x 10- 7 


Br 


3.0 x 10~ 6 



states interaction may not be important. So we think their effects are small and ignore them 
in our calculation. 

Unfortunately, there are no data for these two decays in experimental side up to now. 
We think that the LHC-b experiment can measure these decays in future. The results can 
test this PQCD approach and show some information about new physics. 

The calculated branching ratios in PQCD approach are sensitive to various parameters 
such as the parameters in equations ()16H23p . The uncertainty taken by m has been argued 
in many papers (TTJl EH] , and it is strictly constrained by B — > tt form factor. In Table J21 we 
show the sensitivity of the branching ratios to change of B® and D distribution amplitude 
functions. It is found the uncertainty of the branching ratio in PQCD is mainly due to u b , 
which characterizes the shape of B® meson wave function. 

Considering most of the uncertainty^, we give the branching ratios of these two decays 
with suitable range of u b and a D . Thus we can give our results: 

>0 , n+„-\ q_i_1.2\ w m-7 ( f B s fo \ ( iKTfe ^cs| 



Br(£r -> D + ti-) = (8.3±^) x 10" M — \j_ub_jj±\ — (2g) 

BrffiO -> D-7T+) = (2.9±°"f ) x 10- 6 ( 1 2 ( ^ Ks| 1 2 . (27) 

v s ' y °- 5; V 230 MeV ■ 240 MeV / v - 0412 " - 224 / 

The CP violation information in decay B S (B S ) — > D ±r K T is very complicated. There are 
four kinds of decays 

g = (D+n-\H\B°) oc V^, h = (D+n-\H\B°) oc V cb V* s , 
g=(D-n + \H\B° s ) <xV ub V*, h = (D-n+\H\B s ) <xV; b V US} 



1 Although the uncertainty taken by CKM matrix elements is large, we will not discuss them in this work, 
since they are only an overall factor here for branching ratios. 
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Table 2: The sensitivity of the branching ratios to change of uj, and dp 





(10 7 ) 


(10 6 ) 




td i dO n+ — \ 
DT{B^ — > D^TX ) 


td i d0 n— + \ 
Br(i5g — ► L> 7r^j 


U.4o 


y.o 


O.O 


U.OU 


S Q 

o.o 


o 
z.y 


u.oo 


( .0 


z.o 


an 


T> r ( do . n+/77-— 

ol ^-£3= — ' U 71 J 


ol \D s — > D 71 J 


0.6 


7.6 


2.6 


0.8 


8.3 


2.9 


1.0 


9.1 


3.4 



which determine the decay matrix element of B° s — > D + n~ and D~tt + , and of 5° — > D~tt + 
and D + tt~. They are already shown in the previous Section. There is only one kind of 
contribution for each of the decay mode, thus there is no direct CP violation for them. 
However there is still CP violation induced by mixing, although they are decays with charged 
final states [I]. 

The time-dependent decay rates for B s — ► D ± 7r =F are given by: 

T D± ^(t) = (l±A CP )^^{l + (S D7T ±AS Dn )smAmt 

+ (C Dn ± AC DlT ) cos Amt} , (29) 

and B s -> D ± tt t by 

T D - T (t) = (1±A CP )- {l-[(^±AS D7r )sinAmt 



+(C D7r ± AC D7r )cosAmt]}. (30) 



Utilizing eq. (12*81) . we can get 



\h/g\' 



Cdtt — A C P — 0, ACo-k — ■ i , , |o • 

1 + \h gy 

2 \h/g\ sin7 cos<5 — 2 |/i/gr| sin<5cos7 

^ = i + ' D ^ = i + IW ' 
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(31) 



In deriving the above formulas we have neglected the small weak phase arg(q/p) = p^p? = 
2X 2 7] < 2°, in Wolfenstein parametrization [TT] . 

We can calculate these parameters related to decays B®(B®) — > D ± n T in our PQCD 
approach. Through calculation, we get: 



ACdtt 



-0.56. 



(32) 



The parameters Sd-?t and ASd-k are 7 related. The results are shown in FigureJH If we can 
measure the time-dependent spectrum of the decay rates of and B®, we can extract the 
CKM angle 7 and strong phase 5 in eq. ftBTj) by Figure JH The parameter C D ^ = A C p = is 
from the fact that there is only one kind of contribution to each of the decays. Any deviation 
from zero, will be a signal of new physics contribution. 

.4F ^-—'-\ 








25 50 75 100 125 150 175 



Figure 4: CP violation parameters of B®(B®) 
(solid line) as a function of CKM angle 7. 



D ± n T : ASd-k (dash-dotted line) and Son 



4 Summary 

Recent study shows that PQCD approach works well for charmless B decays [TOl [TT] . as 
well as for channels with one charmed meson in the final states [5j [12] . Because the final 
state mesons are moving very fast, each of them carrying more than 2 GeV energy, there 
is not enough time for them to exchange soft gluons. So we can ignore the soft final states 
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interaction. Due to disadvantages in other approach such as general factorization approach 
[T£] and BBNS approach [19J, pure annihilation decay can be calculated reliably only in 
PQCD approach. 

In this paper, we calculate B® — > Dn decays, which occur purely via annihilation type 
diagrams. The branching ratios are still sizable at the order of 10 -6 . There will also be 
sizable CP violation in these decays. They will be measured in future LHC-b experiment, 
which may bring some information about new physics to us. 



A Some functions 

The definitions of some functions used in the text are presented in this appendix. In the 
numerical analysis we use one loop expression for strong coupling constant, 

An 

a - w = Aw/A 2 )' (33) 

where /3q = (33 — 2nf)/3 and n/ is number of active flavor at appropriate scale. A is QCD 
scale, which we use as 250 MeV at rif = 4. We also use leading logarithms expressions for 
Wilson coefficients presented in ref.jE]. Then, we put m t = 170 GeV, mw = 80.2 GeV, 
and m h = 4.8 GeV. 

The function E\, E m , and E' m including Wilson coefficients are defined as 

E){t) = [d(t) + C 2 (t)3]a s (t) e-^W-s-M (34) 
E m (t) = C 2 (t)a s (t) e -SB(t)-s D (t)SAt)^ (35) 
E' m (t) = Cx(t)a s (t) e-^W-^W-^W, (36) 

where Sb, Sd, and S n result from summing both double logarithms caused by soft gluon 
corrections and single ones due to the renormalization of ultra-violet divergence. The above 
Sb,d,tt are defined as 

S B (t) = 8 ( Xl P+, h) + 2 f ^ 79 (/A (37) 

Ji/bi ^ 

S D (t) = s(x 2 P 2 + , b 3 ) + 2 f ^ 7 ,(A*0. (38) 

Ji/b 2 ^ 

S 7T (t) = s(x 3 P+,b 3 ) + s((l-x 3 )P+,b 3 ) + 2 f (39) 

^1/63 V 
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where s(Q,b), so-called Sudakov factor, is given as 



s(Q,b) 



* Q dn' 

'1/6 

J 67 vr 2 10 2 lB \ faMV Q 



/ ^[{|(2 7E -l-log2)+^log^) 

(40) 



je = 0.57722 • ■ ■ is Euler constant, and 7 g = a s /ir is the quark anomalous dimension. 

The functions h a: ha \ and in the decay amplitudes consist of two parts: one is the 
jet function S t (xi) derived by the threshold resummation [20J, the other is the propagator of 
virtual quark and gluon. They are defined by 

K(x 2 ,x 3 ,b 2 ,b 3 ) = S t (l-x 3 ) (jpj H^{M B ^{l-r^)x 2 {l-x 3 ) b 2 ) 
x {h£\m b ^(1 - r 2 )(l - x 3 ) b 2 )J (M BV / (l - r 2 )(l - x 3 ) b 3 )6(b 2 - b 3 ) + (b 2 <-> b 3j 



(41) 



h{ J \x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = 

|yHf j (MbV(1 - r 2 )x 2 (l - x 3 ) & 1 )J (M B V(1 - r 2 )x 2 (l - x 3 ) b^h - b 

KoiMsF^h), for F 2 } > 



?H«(M B J^ &1 ), for F 2 } 



where Hq 1 ^) = Jo(^) + iY (z), and Fy^s are defined by 

F 2 1} = (1 - r 2 )(x! - x 2 )(l - x 3 ), Ff 2) = xi + x 2 + (1 - r 2 )(l - a* - x 2 )(l - x 3 ). (43) 

We adopt the parametrization for St(x) of the factorizable contributions, 

S t (x) = 2l+ ^^ 2 + +c) [x(l-x)r, c = 0.3, (44) 

which is proposed in ref. [21]. In the non-factorizable annihilation contributions, S t {x) gives 
a very small numerical effect to the amplitude. Therefore, we drop S t (x) in and ■ 
The hard scale t's in the amplitudes are taken as the largest energy scale in the H to kill 
the large logarithmic radiative corrections: 

t\ = max(M BV /(l-r 2 )(l-x 3 ), l/& 2 , 1/63), (45) 



t 2 a = max(M sv /(l-^ 2 )^, I/62, 1/63), (46) 



4 = max(M BA /|F 2 ) |,M sv /(l - r 2 )a; 2 (l - x 3 ), l/&i, I/62). (47) 
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